Pre-classical solutions of the vacuum Bianchi I loop quantum cosmology 
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Loop quantization of diagonalized Bianchi class A models, leads to a partial difference equation as 
the Hamiltonian constraint at the quantum level. A criterion for testing a viable semiclassical limit 
has been formulated in terms of existence of the so-called pre-classical solutions. We demonstrate 
the existence of pre-classical solutions of the quantum equation for the vacuum Bianchi I model. 
All these solutions avoid the classical singularity at vanishing volume. 

PACS numbers: 04.60.Pp,98.80.Jk,98.80.Bp 
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Loop Quantum Gravity (LQG) is the leading candi- 
date for a manifestly background independent approach 
to constructing a quantum theory of gravity y|. This 
approach is particularly well suited in the context where 
the Einstein theory indicates occurrence of singularities 
entailing highly dynamical geometries with extreme cur- 
vatures. The methods employed in this approach can be 
adapted and tested in the simpler context of cosmolog- 
ical models. Quantizing the cosmological models along 
the lines of LQG has lead to the development of Loop 
Quantum Cosmology (LQC) Q. 

One of the crucial simplification available in LQC is 
the existence of the triad representation and knowledge 
of complete spectrum of the volume operator so crucial 
for quantization of the Hamiltonian constraint. The fact 
that the holonomies of the connection are well defined op- 
erators but not the connection it self, is directly respon- 
sible for the two main features of LQC: (a) the quantum 
Hamiltonian constraint leads to a difference equation |3( 
and Op) inverses of triad components have bounded spec- 
tra 0, 0] ■ Both these features lead to the absence of 
'singularities' in the quantum theory 0, Hi • 

While loop quantization of cosmological models is well 
specified, one also needs to check if the quantum dy- 
namics admits solutions (semiclassical states) which can 
approximate the classical description. A natural way to 
recover classical behavior would be in terms of the ex- 
pectation values of suitable observables in the semiclas- 
sical states. LQC (and LQG) being constrained systems, 
makes such a recovery of classical behavior, more compli- 
cated. The solutions of the Hamiltonian constraint, the 
only relevant one for LQC, are typically distributional 
and one needs to equip the space of solutions with a new 
(physical) inner product. One also needs to identify suit- 
able (Dirac) observables. Addressing these aspects is at 
a preliminary stage 0, 0, 0, El . 

Current understanding of the semiclassical limit of 
LQC is centered around the notion of pre-classicality. 
The articulation of this notion has undergone a few 
changes and it is useful to note them. Originally pro- 
posed in [J^ , pre-classicality was thought in terms of ob- 
taining the continuum Wheeler-De Witt equation from 
the fundamental difference equation by a limit in which 
the Barbero-Immirzi parameter is taken to zero. Sub- 



sequently, the idea of pre-classical limit was replaced by 
the idea of a pre-classical approximation To obtain 
(modified) Einstein dynamics from the Hamilton- Jacobi 
equation provided by a WKB approximation [lii H?l Hsj ] . 
it was shown to be sufficient to have approximate pre- 
classical solution(s) with a finite domain of validity. The 
property of local stability is needed in the construction 
of such solutions [l^ and is satisfied by difference equa- 
tions of LQC 0- Recently, the methods of analyzing the 
asymptotic behavior of exact solutions of the difference 
equations have been developed 0, H3| which provide the 
sharpest yet formulation of pre-classicality of a solution 
[lflj . Briefly, the criterion is that asymptotically, small 
scale oscillations in the solution be suppressed. 

Pre-classical solutions so identified are known to exist 
for isotropic models 0,0, H3| and some of the LRS mod- 
els |l9j . However, for the anisotropic, vacuum Bianchi 
I model, pre-classical solutions were shown not to ex- 
ist [2l| . While this result is true for solutions that pass 
through vanishing volume, there are more possibilities 
which permit existence of pre-classical solutions. Since 
the works in 0,01, the so called Bohr quantization has 
been developed [7J which is crucial for the existence of 
pre-classical solutions and we incorporate it in the brief 
summary of the quantum theory given below. 

The kinematical Hilbert space is spanned by orthonor- 
malized vectors labeled as fi2, /is), fii G R. These 
are properly normalized eigenvectors of the triad opera- 
tors p 1 with eigenvalues ^fp/zj, where 7 is the Barbero- 
Immirzi parameter and l\ := SttGH := kH. The volume 
operator is also diagonal in these labels with eigenvalues 
V(/t) given by (^7^p) 3 ^ 2 \/|/ i i/ i 2M3|- Here we have used 
the vector notation to denote the triple (pi, /Z2, ^3)- Im- 
posing the Hamiltonian constraint operator on general 
vectors of the form \s) — X)/i s (aOIm) leads to the fun- 
damental difference equation for the wave function s(fT). 
Here the sum is over countable subsets of R 3 . The wave 
function s(fl) have to be invariant under simultaneous 
reversal of signs of a pair of /xj's and corresponds to the 
residual freedom of reversing the sign of any two of the 
triad components 0. 

In the present context of vacuum Bianchi I model, the 
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Hamiltonian constraint leads to the equation |9j], 

A 12 (p; £12 Mm; £12) + cyclic = , where (1) 

e*i2 = (ei, £2, 4i £2) with each of the e* taking values ±1; 
s(fl;e 13 ) = s(/xi -/i ei -/Xoe'i,M2 - Mo£2 -Mo4iA*3); Mo 
is an order 1 parameter and, 



^4i2(M;£i2) = V(/?;ei2)d(/i3)(eie 2 + 44) 



(2) 



d(/x) 



|1 + a*oM _1 |* - |1 




-1 1 + 



^0 
M = 



Note that while the equation is defined for all 
Mj G K, it is actually a difference equation since only 
the coefficients s(m), differing in steps of ±2mo are con- 
strained by the equation. One can make this explicit by 
setting m/ := 2mo^/ + 2(1qNi, Nj g Z and S^(v) := s(/t). 
Clearly, there are infinitely many 'sectors' labeled by v 
with vj G [0, 1). Only those sectors for which at least one 
of the vi is zero, one will encounter zero volume. 

It is convenient to absorb the factors of volume eigen- 
values into the wave functions by defining t(p) := 
F(m)s(m) so that the equation (JIJ becomes, 

V"<2(M3)(eie2 + 44)*(M;ei 2 ) + cyclic = (3) 



This preserves the gauge invariance condition on the 
s(m). Furthermore, due to the explicit volume eigenval- 
ues, i Ml , M2 . A i 3 — if any of the fii's equal zero. 

Making a product ansatz for the wave function and 
introducing the difference operator A, 



t(p) := Zi(mi) z 2 {[ii) z 3 (m3), (4) 
Az/(m/) := {^/(m/ + 2mo) ~ ^/(M/ - 2 Mo)} , (5) 

the difference equation © can be written as, 

c!(M3)Azi(mi)Az2(m2)z3(M3) + cyclic = . (6) 

The gauge invariance conditions then translate to 
M/) = V z i(j JI 'i)iVI) where 77 = ±1. The vanishing 
of £(/t) when any of the mi = 0, translates into the con- 
dition: z/(0) = 0, VI. 

An exact solution of the partial difference equation can 
be obtained by setting, 

Az/(m/) = fad{Hi)zi{pi) , VI , Vm/ G R (7) 
where /3/ are some constants which have to satisfy, 

Pifo+fofo+Mi = = g>) a -E#- ( 8 ) 

Thus we obtain a class of exact solutions of the partial 
difference equation @ from those of three, ordinary dif- 
ference equations Q with parameters (3i satisfying ©• 



The original partial difference eq. is linear and with 
rea/ coefficients, Atj's. So, without loss of generality, we 
can assume the wave function, s(/Z) to be real. Since for 
the product ansatz, all zi are independent, these must 
be real as well which requires that the /3/'s be real. 

If all (3j are zero, then zj's are constants and so is t(fx). 
These cannot satisfy the condition zj(0) = without 
making t(ft) — identically. The non-trivial solutions 
then cannot pass through zero volume and must belong 
to the sectors with vj ^ 0. This also shows that in the 
sectors with v\ ^ 0,VI, there is an exact pre- classical 
solution to namely t(2[iQV + 2/ioN) — constant 
0),V7V G Z 3 . 

If at least one /3j ^ 0, then we can always take out 
a common factor from all /3j and ensure Yli 01 = 1- 
Equivalently, a common scaling of /3j can be absorbed 
by a common inverse scaling of the d(p) functions, which 
amounts to a scaling of the volume which cancels out in 
The class of solutions that is being constructed can 
thus be parameterized exactly in the same manner as the 
classical Kasner solution. In particular either two of the 
P's are zero or exactly one is negative while other two are 
positive. From now on we will restrict to < \(3i\ < 1. 

To explore pre-classicality of the separable solution, 
let us focus on (J2J), suppressing the label I. Introducing 
the notation: m ;= 2mo^ + 2mo", n G Z, v G [0, 1) and 
z(m) := Z„ , can be written as: 



K+2 = Z v n + pd{u,n+l)Z v n+1 
1 



d(u, n) 



1 



2{v 



1 



Vn G Z (9) 

1 

1 



2(v + n) 



We have infinitely many decoupled sectors, labeled by 
v and for each of these we have a second order, ordinary 
difference equation. Due to linearity, only one condition 
is enough to determine a solution. Only for the sector v = 
(m is integer multiple of 2mo), the condition z(0) = is 
relevant and it fixes the solution completely. The gauge 
invariance condition translates into the identification: 



(n+i) = vZn, V n G Z, 77 



±1 



(10) 



which restricts the Z's in the same sector, only for 
v = 0, |. For all other sectors the gauge invariance condi- 
tion relates two different sectors. Under the identification 
implied by IjlOp. the equation satisfied by Z^ goes over 
to the equation satisfied by Z]^ v automatically with the 
same value of 0. Therefore the solutions in the sector 
(1 — v) can all be obtained from solutions in the sector v 

via JTDJl. For v=\, the gauge invariance condition re- 

1/2 1/2 

quires Z_ l = r/Z which fixes the solution completely. 
For the two sectors, v = 0, |, pre-classicality is not op- 
tional - the solution is either pre-classical or it is not. 
From now on the superscript v is suppressed. 

Consider Defining the ratios u n +i ■— Z n+ i/Z n , 
one can see that [(u n +2 — (3d{v,n + — l]Z n = 0. 
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There are several possibilities now (n > for definite- 
ness). If Z n remains non-zero for all n, then u n — > ±1 as 
n — » oo. The u n — > — 1 is referred to as a sequence with 
sign oscillations. Since u n+ 2 — /3d(u,n + 1) and u n+ \, 
must have the same sign, it is clear that if Ui > and 
f3 > 0, then u n > for all n. For /3 < however, u n may 
become negative for some no and then stay negative sub- 
sequently. Whether this could happen depends on the 
value of u\. Thus it is conceivable that for some positive 
values of u±, one could have a sequence without sign os- 
cillations. If Z n converges to zero, then \u n \ converges 
to a value < 1, once again allowing for u n — » — 1. The 
sequence u n could also converge to in which case there 
could be oscillations about 0, but these are suppressed. 
By pre-classicality we mean either absence or suppression 
of sign oscillations. 

To identify pre-classical sequences, we employ the gen- 
crating function technique [lflj . The function being 
an algebraic function poses some difficulties which can 
be handled in exactly the same way as in [Tjj ■ Basically, 
one separates out the large n part of d(i>, n) and solves 
the equation perturbatively, Z n = a n + Yl'uLi %n- The 
leading order term, a n satisfies the equation with d(v, n) 
replaced by (2(n + v))~ l - Let us focus on n > so that 
the absolute signs can be removed. Setting = 4A the 
equation defining the a n sequence is, 



2A 



a n +2 



1 







n > 



(11) 



To account for non-integral v, we define a generating 
function F(x) := J2^=o a nX n+u and the function G(x) :— 
x~ 1 (F(x) — aox u ), which satisfies a differential equation 
equivalent to the difference equation Ijllfl. 

-^-((l-x 2 )G)-2\G-a (v+l)x u -axvx v ~ x = (12) 
ax 

For non-zero u, the last term is singular at x — im- 
plying that G{x) will not be analytic at x = 0. However, 
the singular term is integrable such that G(x) is contin- 
uous at x = 0. The series representation requires it to 
vanish at x = 0. 

The equation l|12l) can be easily integrated to give, 



G(x) 



fl 



-x-i 



(13) 



Notice that given any <Zo,ai, the solution to (|ll|l is 
completely determined and so should be G(x). The above 
solution for G(x) contains an indefinite integral and an 
arbitrary constant of integration, cq. We must choose cq 
and convert the indefinite integral to a definite one such 
that G(x) corresponds to the sequence specified by the 
given ao,ai. The only value of G(x) we know without 
having to know the full sequence is G(0) = 0. Further- 
more both the integrands in the integrals in (|13J) are in- 
tegrable at x — 0. Thus it is possible to impose G(0) = 



which then requires Cq 
then obtained as, 



0. The generating function is 



(1 - x)G(x) 
l(v, A, x 



(l + a;) A - 1 (l-x)- A x (14) 
la I(v, A, x) + axl{v - 1, A, x)] , 
1-f 



) := + 



1 + t 



i" 



(15) 



One may already note that the integrals are finite at both 
x = ±1 but the pre- factor is not. Thus singularities 
of (1 — x)G(x) are controlled by the pre-factor which 
is independent of v. Apart from the v dependence of 
the integrals, one does not expect qualitative behavior of 
G(x) to be affected by v. Furthermore, since there are 
two free parameters (ao,di) which specify the sequence 
and only one of these is relevant one due to the linearity, 
one can at the most impose only one condition capturing 
pre-classicality, to get a non-trivial solution. 

Now consider the behavior of G(x) as x — > —1. For 
the range of < |A| < 1, both the integrals exist (and 
are positive), but the pre-factor diverges. A divergence 
in G{x) at x = —1 implies un-suppressed sign oscillations 
which are to be avoided for pre-classical sequences [19|. 
Clearly, to avoid this singularity in G(x), the integrals 
must add up to zero which determines a\ in terms of do- 
For x = — 1, one has, 

i(u,\,-i) = {-\yi{v,-\i) , (i6) 

I(v, -A, 1) = (1 + v)B(l + v, 1 - A) x 

F(-X, v + 1, -A + v + 2; -1) , (17) 

where, B and F are the Beta function and the hypergeo- 
metric functions |22l| . The condition of no singularity at 
x = — 1 gives, 



aoI(v, —A, 1) = a\l(v — 1, —A, 1) 



(18) 



This determines the sequence uniquely modulo a triv- 
ial, overall scaling. The sequence satisfying (|18fl has sup- 
pressed sign oscillations. Its convergence properties are 
determined by the x — ► 1 behavior of (1 — x)G[x). For 
x = 1 the integrals again exist but now the pre-factor 
diverges for A > and vanishes for A < and so does 
the sequence {a n }. However, (1 — x)G(x) is integrable 
at x = 1 which implies that the asymptotic behavior 
of a n is bounded by n. Making a power law ansatz for 
asymptotic a n , one can see from (jl 1|> (and indeed from 
© as well) that a n ~ n x . Since (3i come with both 
signs, both behaviors must be admissible. For A = 0, one 
gets a\ = ao and the sequence is the constant sequence 
a« = cio Vrt > which is obviously pre-classical. 

As anticipated, these results are exactly analogous to 
those obtained in 0|. Indeed, for v = 0, equations (|lll 
IT21 Ibljl go over to the equations of 0] . Now the bound- 
ary condition is G(0) = a\ which gives cq = a\ and def- 
inite integral has the lower limit as 0. Demanding non- 
singularity of G(— 1), determines a\ exactly as in |l9| . 
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The behavior at x = 1 is similar to that for the non-zero 
v case. 

Thus, there certainly exist sectors such that for each 
choice of the separation constant f3, < |/3| < 1, one 
can select a unique solution of which is pre-classical. 
These solutions of course have to be improved by com- 
puting the corrections Z^ [lflj ]. The asymptotic power 
law behavior of a n will continue to hold also for Z n . 

All these statements hold for n > 0. Having deter- 
mined Zq,Z", by pre-classicality, Z n< o can be deter- 
mined by the exact equation 0. Whether Z% is pre- 
classical also for n < 0, can be inferred by testing for 
pre-classicality of Z]^~ u for positive n, using (fTrj|) . 

As noted earlier, the sectors v = 0, \ already have a 

unique solution due to the conditions Zq = and Z]_^ — 

1/2 

i]Z respectively. If these conditions are imposed on the 
leading pre-classical sequence {a n }, then clearly there are 
no non-trivial solutions in the v = sector. For v = 1/2 
sector, numerically, the gauge invariance condition and 
the pre-classicality condition seem to hold only for (3 = 
with rj = 1. For other sectors, pre-classicality is the only 
condition imposed and solutions can be constructed. 

The full wave function is the product of the three se- 
quences and apart from an over all constant factor, is 
completely determined. One can build more general (and 
non-separable) solutions by taking complex linear com- 
binations with coefficients being functions of (3. Clearly, 
a combination involving a diverging and a vanishing so- 
lution will be a diverging one and still without sign oscil- 
lations. Since such a solution involves a distribution of 
j3, these parameters themselves would not be identified 
with the classical Kasner parameters, ai (say) satisfy- 
ing J2i a i = 1 = J2i a i- Rather, one would imagine 
constructing linear combinations which 'peak' in some 
suitable sense, around three triad values pg and a Kasner 
parameter q?o- (Since the reduced phase space of the vac- 
uum Bianchi I model is 4 = 6 — 2 dimensional, one needs 
four parameters to specify a classical solution and these 
could be conveniently taken as three initial triad values 
and a Kasner parameter value.) If such a construction 
can be carried out, then one would be able to claim that 
the quantum theory has 'sufficient number of scmiclassi- 
cal solutions' as expected from the classical theory. 

We note that since solutions of the Hamiltonian con- 
straint are expected to be distributional in general, kine- 
matical normalizability of the pre-classical (or otherwise) 
solutions is not directly mandated. The requirement of 
pre-classicality for both signs of n is an open issue. For 
an alternative treatment of separable solutions, see [24| ■ 

In summary, we have shown that in every sector vj ^ 
0, ^, there exist a one parameter family of pre-classical 
solutions. For fii = 0, the solution is in fact exact, possi- 
bly corresponding to the Minkowski space-time. All these 
solutions skip the vanishing volume eigenvalues. By con- 
trast, in the vi = sector, there are no pre-classical so- 
lutions [2ll | . The richness of the loop quantization, man- 



ifested by the infinitely many sectorSjis crucial for this 
result; an observation also made in hll in the isotropic 
context. We have heuristically indicated how these fam- 
ilies can be used to see if loop quantization does admit 
'enough semiclassical states'. The exact, non-singular 
solutions of the effective dynamics of vacuum Bianchi 
I models given in j2^| also exhibit a similar feature of 
avoiding vanishing volume which motivated this work. 
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Hossain are also acknowledged. 
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